Effect of Long-Range Interactions in the Conserved Kardar-Parisi-Zhang Equation 
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The conserved Kardar-Parisi-Zhang equation in the presence of long-range nonlinear interactions 
is studied by the dynamic renormalization group method. The long-range effect produces new 
fixed points with continuously varying exponents and gives distinct phase transitions, depending on 
both the long-range interaction strength and the substrate dimension d. The long-range interaction 
makes the surface width less rough than that of the short-range interaction. In particular, the 
surface becomes a smooth one with a negative roughness exponent at the physical dimension d = 2 



PACS No.: 05.40.+j; 05.70.Ln; 68.35.Fx; 6I.50.Cj. 

For the last decade the kinetic roughening of surfaces 
has attracted much interest B. The recent studies con- 
centrate on measuring the scaling exponents which char- 
acterize the asymptotic behavior of the surface roughness 
on a large length scale and in a long time limit, and find- 
ing the continuum equations. The problem of rough sur- 
face is not only of practical importance in crystal growth, 
but also related to the nonequilibrium statistical physics. 
Lots of computer simulations and theoretical approaches 
have been applied in the studies of the Kardar-Parisi- 
Zhang (KPZ) Q equation and discrete molecular-beam- 
epitaxy (MBE) growth models with various kind of noises 
. Among them, the Eden model J8| , ballistic deposi- 
tion ||, and the restricted solid-on-solid (RSOS) growth 
model [|l0| have been identified as a universality class cor- 
responding to the KPZ equation for the coarse grained 
height variable h(r,t) which describes the surface as a 
function of coordinate r and time t. The KPZ equation 
has a nonlinear term of short range describing the lateral 
growth. However, there is a poor agreement between the 
KPZ equation and the experimental data. 

Recently, Mukherji and Bhattacharjee proposed a 
phenomenological equation in the presence of long-range 
interactions to describe the kinetic roughening of the sur- 
face growth. The long-range effect of the nonlinear term 
in the KPZ equation is introduced by coupling the gradi- 
ents at two different points. The roughness of the surface 
are found to depend on the long-range nature and several 
distinct phase transitions are observed. The long-range 
interactions decaying slower than 1 /r d (d is the substrate 
dimension) makes the KPZ fixed point with the short- 
range interaction be unstable. The surface then have the 
long-range roughness with different exponents depending 
on the power law of the long-range interactions. Other 
interactions decaying faster than 1 jr d are suppressed by 
the local interaction yielding the ordinary KPZ univer- 
sality class. 

In the kinetic roughening problems, the universality 
class of the dynamic systems depends on the symmetry 
of the order parameter, the dimensionality of space, and 



the conservation of the surface currents. Therefore it 
would be interesting to examine how the long-range in- 
teraction in the conserved growth equation affects both 
the roughness of the surface and the phase transitions 
compared to the cases of the short-range interaction. We 
extend the phenomenological equation of Mukherji and 
Bhattacharjee, to a conserved equation, 



dh(r, t) 
dt 



= -KV 4 h(r,t)+ Vc (r, t) 



-V 2 J dr'-d(r')Vh{r + r',t) ■ Vh{r-r',t) , (1) 

where h(r,t), assumed to be a single- valued function 
of position r, describes the height of the surface. The 
parameter if is a constant, and r\ c is a conserved 
random noise of zero mean with {r) c (r 7 t)rj c (r' , t')} = 
-2D c V 2 S{r - r')S(t - t'). Since the right hand side of 
Eq. (|l]) can be written as the divergence of a current, 
the total volume under the surface is conserved. The 
kernel ?9(r) includes long-range part which is connected 
to the underlying interactions. As Ref. O], we take 
$(r) to have a short-range(SR) part Ao<5(r) and a long- 
range(LR) part ~ r p ~ d , or more precisely, in Fourier 
space, t9(k) = Ao + \ p kr p . 

The surface width W(L, t) can be described by the dy- 
namical scaling form, W(L, t) = L x F(t/L z ), where L, x, 
z and F are the system size, the roughness exponent, the 
dynamic exponent, and the scaling function respectively. 
For Ao = A p = 0, it becomes a linear equation evolving 
with the conservative surface diffusion, where the rough- 
ness exponent x is (2 — d)/2, and the dynamic exponent 
z is four. For the physical dimension d = 2, x = 5 
thus the surface width is logarithmically rough as a func- 
tion of system size L. Above two dimensions the lin- 
ear equation with the conservative noise produces nega- 
tive roughness exponent, implying a smooth surface |JL2j. 
For A p = and A ^ 0, Eq. (|l|) becomes the conserved 
KPZ equation with a conservative noise (called Sun-Guo- 
Grant(SGG) equation) [Q,^3|, where the average height 
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remains constant. For this local conserved growth equa- 
tion, the dynamic renormalization groupfRG) calculation 
shows x = (2-d)/3 and z = (10 + d)/3 fj. For d > 2, the 
nonlinear term is irrelevant and then the exponents are 
given by the linear theory with both Ao and A p being zero 
in Eq. (Q). Here we show that a long-range part(A p ^ 0) 
gives a new fixed point with continuously varying expo- 
nents and thus yields distinct phase transitions depend- 
ing on both the parameter/? of the long-range interac- 
tions and the substrate dimension d. This nonlocal A p 
term with positive p makes the surface less rough than 
the case of A p = 0. Especially, at the physical dimen- 
sion d = 2, the surface becomes a smooth phase with a 
negative roughness exponent rather than logarithmically 
rough phase as in the SGG case. 

Under the change of scale, the parameters in Eq.([l]) 
change to K -> b z ~ 4 K, D r -> b z ~ 2 *- d - 2 D r , A n -> 



and R = Uq/U p , we obtain the flow equations for Uq, U p , 
and R: 



& z+ x- 4 A , and A, 



y-+x+p 4\ j n ^ ne aDsence f 



nonlinearity(Ao = X p = 0), K and D c are scale invariant 
to yield z = 4 and xo = (2 — d)/2. Using these values we 
find that the nonlinearities rescales as Aq — > 6' 2_d ^' 2 Ao 



and A p -» b<- 2 + 2 P- 



'\ p . So that the critical dimensions 



are given by d c = 2 + 2p (p > 0) and d c — 2 (p < 0) for 
any nonzero A p . When p > 0, if d < d c = 2 + 2p, the 
fixed point of the local interaction (A p — 0, Ao ^ and 
2 + x ~ 4 = 0) is unstable and thus a new fixed point is 
expected. If 2 + 2p < d, the nonlinearities become irrele- 
vant and the surface is controlled by the linear equation. 
For p < 0, if d < 2, the SGG fixed point is stable so that 
Ao is relevant rather than A p , otherwise the linear term is 
relevant. As a result, various phase diagrams depending 
on d and p would be appeared. 

Following a dynamic RG procedure integrating 
out fast modes in the momentum shell e A < |k| < A 
and performing the rescaling r — > br, t — > b z t, h — > b x h, 
we derive the following flow equations for the coefficients, 
in a one-loop approximation, 



dK 

~dt 

dPg 

dt 
dX 

dt 
dX p 

dt 



K 3 

D c [z-2 X -d-2] 
A [z + x-4] 
A p [z + x~4 + p] 



Ad 



(2) 
(3) 
(4) 
(5) 



where f{a) = d\n§{k) / d\nk\ k=a , and B d = S d /{2^) d , S d 
being the surface area of a d-dimensional unit sphere. 
Since the diagrams contributing to D c have prefactors 
proportional to fc 4 , they correspond to higher derivatives 
in the original noise spectrum. Note that two scaling re- 
lations, z + x = 4 and z + x — 4 — p which result from 
the non-renormalization of the Ao and A p in Eq. (|4j) and 
Eq.(||) respectively, are the results of one-loop approxi- 
mation |15| . 

Defining the dimensionless 

parameters U 2 = {D c X 2 B d )/K 3 , U 2 = (D c X 2 p B d )/K 3 , 



dUo 
dt 
dU, 



= U 



2-d , 3{d-A) Tr2 t 3C/ P 



Sd 



U^—ficoUo+cxUp) 



—?- = U 

dt p 



2-d + 2p t 3(d-4) 2 



(6) 



(7) 



and dR/dt = -pR, where Co = (d - A)2- p + d - A - 3p, 
and ci = (d — 4 — 3p)2 _p . The equation for R rules 
out the existence of any off-axis fixed point in the Uq 
and U p parameter space (except for p = 0). From these 
equations we find that there are only two sets of axial 
fixed points in the two dimensional (Uq, U p ) space: SR = 
{C/ * 2 = Ad(d-2)/3(d-A),U; 2 = 0}, with x + z = 4, and 
LR = {US 2 = Q,U; 2 = Ad(d-2-2p)/3(d-A~3p)2-P}, 
with x + z — 4 — p. When U p = 0, the SR fixed point is 
stable for d < 2, where X = (2 - rf)/3 and z = (10 + d)/3, 
in agreement with the results of Sun, Guo, and Grant [Q. 
For d > 2, C/o is driven to zero as I — > oo. The surface 
width is thus described by the linear equation yielding 
a smooth phase except for d = 2 (logarithmically rough 
phase). Similarly, from Eq. (0) with Uq = 0, the LR 
fixed point for d < 2 + 2p is stable. At this new LR fixed 
point, the exponents are given by 

X=(2-d-p)/3, and z = {10 + d - 2p)/3. (8) 

These exponents are determined by Eq.(|^) and Eq.(j^) 
in which D c and A p are not renormalized in an one-loop 
approximation (z — 2x — d — 2 = and z + x = A — p, 
respectively). 

From these recursion relations, we can discuss the sur- 
face morphologies and the phase transitions for all d's 
and p's (see FIG. 1). Note that Eq.(]l|) is invariant under 
h — > —h and A — > —A transform. Therefore we consider 
both positive and negative values of U p and take Uq > 
without any loss in generality. As shown in FIG. 1, there 
are various (U p , Uq) phase diagrams depending on the di- 
mensionality d and the long-range interaction parameter 
p. We explain the various phase diagrams in detail. 

I) p > 0: The effective nonlinearity U p is dominant 
than C/o, then the phase in all space (U p ,Uq) except 
for U„ — is determined by the long-range A p term in 
Eq. (|l|). For d < 2 — p, the LR fixed point is stable, and 
the surface is the LR rough phase with the positive rough- 
ness exponent (we call it the LR rough) given by Eq. (ra). 
If Up = 0, the SR rough phase with the positive roughness 
exponent (x = (2 — d)/3, we call it the SR rough) exists 
such that a phase transition takes place between two LR 
rough phases when the sign of U p is changed. The critical 
behavior(£/ p = 0) follows the SGG's nonlinear equation 
with the SR rough phase. For d — 2 — p, the surface is the 
logarithmically rough phase with a zero roughness expo- 
nent (we call it the Log rough). For 2 — p < d < 2 + 2p, 
the phase is controlled by the LR fixed point and the 
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FIG. 1. U p (y-axis) vs (7 (x-axis) phase diagram in (p, d) 
space. On the a;-axis and y-axis of the phase diagrams, the 
solid lines denote the rough phase and the dotted lines do a 
smooth phase. The detailed meanings of the corresponding 
rough or smooth phase are explained at Table 1. For p > 0, 
the long-range effect makes the surface less rough than a case 
of p — (see, the region 2 — p < d < 2 for p > 0). 

the surface is the LR smooth phase due to the negative 
value of the roughness exponent (x = (2 — rf — p)/3, we 
call it the LR smooth). In addition, the various critical 
behaviors depending on both value of p and dimension d 
are shown in FIG. 1. Phases at the critical \ine(U p = 0) 
are SR rough for 2 — p < d < 2, Log rough for d — 2, and 
MH smooth for 2 < d < 2 + 2p. Here, the MH smooth 
phase is defined by the linear equation with a negative 
roughness exponent (\ = (2 — d)/2, we call it the Mullins- 
Herring(MH) smooth). For d > 2 + 2p, both the LR 
and SR fixed points are irrelevant, so that only the MH 
smooth phase of the linear equation exists. Therefore, 
various phase transitions take place when the sign of U p 
is changed, except for the region (d > 2 + 2p) where 
no phase transition occurs for all values of U p and l7 . 
At physical dimension d = 2, it is well known that, for 
the short-range interaction)^ = and Uq ^ 0), the 
SR fixed point(SGG) is irrelevant, so that the surface is 
logarithmically rough. However, if U p ^ and p > 
(that is, for the long-range interaction), the LR fixed 
point is relevant and the surface becomes LR smooth 
with the negative exponents given by Eq.(^). We thus 
find that nonzero U p term with p > can make the 
surface less rough than the logarithmically rough of the 
case p = (see Table 1). 

II) p < : The LR fixed point is irrelevant on the 
ground that Un is dominant than U p . So, the short-range 
term in Eq. (yj) corresponding to the SGG equation de- 
termines the surface behavior in all space (U p , Uo) except 
for U = 0. For d < 2 + 2p, the SR fixed point is stable 
and the surface is always SR rough except for Uo = 0, 
when it is a LR rough phase. So, there is no phase tran- 
sition for U ^ 0. For d > 2 + 2p, both the SR and LR 
fixed points are no longer stable, so that the phase is 
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Log rough 
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MH smooth 


MH smooth 
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SR rough 


LR rough 
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SR rough 


Log rough 




2 - p < d < 2 


SR rough 




p > 


d = 2 


Log rough 


LR smooth 




2 < d < 2 + 2p 


MH smooth 






d> 2 + 2p 


MH smooth 


MH smooth 



Table 1. Various phases depend on both p and d. These 
phases correspond to the diagrams in FIG. 1. There are six 
different phases: LR rough if \ — (2 — d — p)/3 is positive, 
LR smooth if \ — (2 — d — p)/3 is negative, SR rough if 
X = (2 — d)/3 is positive, MH rough if % = (2 — d)/2 is posi- 
tive, MH smooth if \ = (2 — d)/2 is negative, and Log rough 
if X = o. 



governed by the linear equation. Therefore, the phases 
are MH rough for 2 + 2p < d < 2, logarithmically rough 
for d — 2, and MH smooth for d > 2. Here, the MH rough 
phase means the positive roughness exponent given by 
the linear equation (x = (2 — d)/2). Unlike the case of 
p > 0, any phase transitions do not take place for all 
spaces of U p and Uo- 

We have also studied Eq.(]l|) with a nonconservative 
noise r\ instead of a conservative noise rj c . The non- 
conservative noise r] is a white noise of zero mean with 
{ri(r,t)rj(r',t')) = 2DS(r - r')S(t - t'). There are two 
sets of axial fixed points. The phase diagrams in this 
growing surfaces are essentially the same as the FIG. 1 
if the dimensionality is replaced by d — > d — 2. At the 
SR fixed point, d c = 4, x = (4- d)/3 and z = (8 + d)/3, 
in agreement with those of the equation introduced by 
Lai and Das Sarma At the new LR fixed point with 
d c = 4 + 2p, the exponents are given by x = (4 — d — p)/3 
and z = (8 + d — 2p)/3 which are obtained from the 
non-renormalization of D and A p in one-loop approxima- 
tion. Experimental results for the growth of Fe films on 
Fe(001) using high-resolution low-energy electron diffrac- 
tion(HRLEED) technique fl6| show x = 0.79 ± 0.05 and 
= x/z = 0.22 ± 0.02 for d = 2. These values are con- 
sistent with our results with the value p — —0.37. At 
this point, it is unclear whether this experimental sys- 
tem really possesses the long-range interaction as that 
of Eq. ([!]). It is thus strongly invited to examine other 
systems with the long-range interactions. 

In summary, we have studied the conserved KPZ equa- 
tion in the presence of long-range interactions. For posi- 
tive value of p, the long-range nonlinear term makes the 
surface less rough and produces the different values of 
the exponents compared to them of the SGG equation. 
Particularly, at physical dimension d = 2, the surface is 
smooth for p > 0, while the surface is logarithmically 
rough for p = 0. However, the long-range nonlinear term 
becomes irrelevant for the negative value of p, and the 
surface is controlled by the SGG fixed points. 
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